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We  have  conducted  research  on  the  modeling  and  control  of  nonlinear  systems.  Our  efforts  have 
been  directed  toward  understanding  the  control  of  truly  nonlinear  behavior  as  well  as  the  synthesis 
of  rantrol  laws  for  systems  that  can  be  nearly  transformed  into  linear  systems  (via  approximate 
feedback  linearisation). 

We  have  introduced  a  new  framework  for  understanding  and  analyzing  the  stability  and  control 
nonlinw  maneuvering  systems.  This  approach  is  based  on  the  concept  of  innsverse  dynamics. 
We  have  demonstrated  the  usefulness  of  this  approach  b  the  control  of  the  swingbg  energy  of 
the  pendulum  for  an  experimental  cart  and  pendulum  system.  On  the  theoretical  side,  we  have 
provided  a  new  method  for  the  construction  of  converse  Lyapunov  functions  for  exponentially 
stable  periodic  orbits.  ' 

New  teclmques  for  the  approximate  feedback  Ibearizatimi  of  nonlbear  systems  have  been 
devdo^.  In  order  to  construct  a  feedback  linearising  coordbate  transformation,  a  class  of 
(^muz^ion  probians  has  been  formulated  for  finding  approximate  sdutkms  to  an  wprtmriate 
systOT  m  partial  differential  equations,  h  contrast  to  prevbus  results,  this  uproaeh  does  not 
req^  dtfferenitaiton  of  the  data  describbg  the  system  and  is  therefore  appUcable  to  systems 
with,  for  example,  tabular  data. 
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We  have  conducted  research  on  the  modeling  and  control  of  nonlinear  systems.  Our  efforts  have  been  di¬ 
rected  toward  understanding  the  control  of  truly  nonlinear  behavior  as  well  as  the  synthesis  of  control  laws 
for  systems  that  can  be  nearly  transformed  into  linear  systems  (via  apprmdmate  feedback  Unearuation). 

Numerical  Approximate  Feedback  Linearization 

We  have  studied  the  approximate  feedback  linearization  of  nonlinear  systems  of  the  form 

X  =  /(*)  +  g{x)u 

where  f  and  g  are  system  vector  fields  that  may  be  available  only  as  numerical  data.  The  problem  is  to  find 
coordinate  functions  <  =  1,  ■ . . ,  n,  so  that  the  transformed  system  =  d>i(x)) 

4i  =  ^2  +  V’i(»)  +  ^i(»)« 

^n-l  =  -I- V>n-l(*)  +  ffn-l(*)« 

in  =  d(x)4-a(a!)«  -|- V’i»(*)  +  ^n(*)« 

is  close  to  a  linear  system  over  some  (compact)  region  D  of  the  state  space.  That  is,  we  wish  to  make 
V'j(x)  ;=  Ljtt>i{x)  —  0f+i(x)  and  ^,(x)  ;=  Lf<t>i{x)  small  over  Q. 

We  have  developed  and  implemented  algorithms  [2,  3]  to  construct  the  necessary  coordinate  functions  by 
solving  an  unconstrained  weighted  least  squares  problem  of  the  form 

r  f"“* 

+  *"#(*) 

where  .F  is  a  suitable  finite  dimensional  function  space.  In  particular,  we  have  selected  a  space  of  tensor 
product  B-splines  defined  over  D  with  useful  boundary  conditions  (e.g.,  (x)  =  xi  and  ^,(x)  =  0,  t  =  2, . . . ,  n 

on  the  equilibrium  manifold  €).  Noting  that  ^i(x)  =  Bj(x)a)  with  basis  functions  Bj{-)  this  problem 
then  becomes 

min  a^Qa  +  p'^a  +  r 
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which  can  be  solved  by  solving  the  linear  system 

2Qa  =  —p. 

The  dimension  of  the  coefficient  vector  a  can  become  somewhat  large  depending  on  the  dimension  of  the 
space  as  we' 1  as  the  dimension  of  the  function  space.  Fortunately,  due  to  the  finite  support  property  of  B- 
splines,  the  matrix  Q  is  sparse.  The  elements  of  Q  and  p  are  computed  by  estimating  the  above  integral  over 
the  small  cubes  associated  to  each  basis  function.  By  fitting  the  components  of  /  and  g  with  polynomials 
over  each  of  these  cubes,  we  have  developed  an  efficient  quadrature-like  technique  for  estimating  the  elements 
of  Q  and  p. 

We  have  seen  that  this  approach  is  effective  in  constructing  feedback  linearizable  nonlinear  system  ap¬ 
proximations  to  a  given  nonlinear  system.  Unlike  most  traditional  feedback  linearization  approaches,  this 
new  technique  does  not  require  the  differentiation  of  system  vector  fields  and  is  thus  applicable  to  systems 
where  the  vector  fields  are  only  available  numerically.  Such  approximations  can  be  used  to  develop  nonlinear 
control  systems. 

Several  important  questions  remain.  For  example,  what  is  a  fair  way  to  compare  the  performance  of 
various  nonlinear  control  schemes?  In  many  cases,  it  is  difficult  to  characterize  what  behavior  we  would  like 
(or  should  expect)  the  nonlinear  system  to  exhibit.  This  is  in  marked  constrast  to  current  state  of  affairs  for 
linear  systems. 

Unfortunately,  most  control  objectives  continue  to  be  specified  from  a  linear  point  of  view.  This,  in  many 
cases,  leads  to  quite  acceptable  results  but  may  be  unnecessarily  limiting  our  possibilities.  In  order  to  break 
free  from  this  line  of  attack,  we  feel  that  it  is  crucial  that  a  number  of  nonlinear  benchmark  problems  be 
developed.  Furthermore,  it  is  essential  that  these  benchmark  problems  involve  a  significant  experimental 
component.  Without  an  experimental  component,  many  important  nonlinear  effects  will  not  be  uncovered 
(we  will  be  merely  stretching  our  linear  ideas). 

Nonlinear  Control  about  a  Periodic  Orbit 

Much  of  the  research  in  nonlinear  control  has  focused  on  extending  linear  system  theory  and  results  to 
nonlinear  systems.  However,  since  linear  systems  are  a  strict  subset  of  class  of  nonlinear  systems,  this 
approach  will  certainly  miss  many  important  possibilities.  With  this  in  mind,  we’ve  been  investigating  the 
control  of  nonlinear  systems  around  period  orbits — the  simplest  truly  nonlinear  phenomenon.  Furthermore, 
to  reduce  the  prevalence  of  artificial  assumptions,  we  have  included  a  substantial  experimental  component 
in  our  investigation. 


Figure  1:  A  cart  and  pendulum  system 

Our  research  has  been  motivated  by  the  common  cart  and  pendulum  system  shown  in  Figure  1.  Rather 
than  the  usual  objective  of  balancing  the  inverted  pendulum,  we  chose  to  make  regulation  of  the  swinging 
energy  of  the  pendulum  the  control  objective.  The  results  of  this  research  were  presented  at  the  Slst  CDC 
and  have  been  accepted  for  publication  [1]. 

The  (normalized)  dynamics  of  this  system  have  the  simple  form 


$  =  -nn0  — coe0u 

ic  =  n 


(1) 
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where  the  input  u  is  taken  as  the  cart  acceleration.  The  objective  was  to  regulate  the  swing  energy  H{0,ui)  := 
w*/2+(l-costf)  (where  w  =  9)  to  a  desired  swing  energy  H.  Defining  the  energy  error  E{0,u)  :=  H(6,u)-H 
we  chose  the  control  law 

U  =  aWCOstfE  —  (2CWeVe  +  (2) 

where  a,  ^d  Uc  design  parameters  and  Ve  =  Xe-  This  control  law  (structure)  was  selected  on  the 
following  grounds.  If  the  cart  dynamics  is  completely  ignored  (in  the  control  and  the  dynamics),  the  closed 
loop  error  dynamics 

E  =  —au)^  cos*  6E 

can  be  shown,  using  the  theory  of  Poincare  maps,  to  be  exponentially  stable  (for  almost  all  initial  conditions). 
The  remaining  portion  of  the  control  law  was  added  to  stabilize  the  cart  position  without  destroying  the 
stability  of  the  periodic  orbit. 

The  stability  and  robustness  characteristics  of  the  closed  loop  system  are  best  understood  by  looking  at 
the  dynamics  transverst  to  the  desired  periodic  orbit.  For  this  system,  E,  Xt,  Vc  provide  suitable  transverse 
coordinates.  By  showing  that  the  (periodic)  time  varying  linearization  of  the  transverse  dynamics  about  the 
periodic  orbit  was  asymptotically  stable,  we  were  able  to  conclude  that  the  orbit  was  itself  exponentially 
stable.  Furthermore,  since  we  analyzed  the  (time-varying)  linear  system  using  a  small  gain  approach,  the 
overall  system  possessed  additional  robustness  properties.  These  characteristics  were  further  confirmed 
experimentally.  The  experimental  component  also  led  to  interesting  insights  concerning  the  use  of  state 
estimators  for  nonlinear  control  [1]. 

In  general,  the  local  dynamics  about  a  periodic  orbit  can  be  described  by 

9  =  l  +  /i(tf./>) 

p  =  A{$)p-Vh{9,p)  ^  ^ 

where  /i  and  are  first  and  second  order,  respectively,  in  the  transverse  coordinates  p.  Many  powerful 
analysis  and  synthesis  tools  can  be  contracted  by  noting  that  the  periodic  orbit  is  exponentially  stable  if 
and  only  if  the  transverse  linearization 

|  =  (4) 

is  asymptotically  stable. 


Figure  2:  A  domain  of  attraction  estimate  for  a  periodic  orbit. 

We  have,  for  example,  used  this  fact  to  provide  a  means  for  constructing  converse  Lyapunov  functions 
for  exponentially  stable  periodic  orbits  [4].  Such  functions  can  be  used  to  estimate  the  domain  of  attraction 
of  the  stable  periodic  orbit.  Figure  2  illustrates  the  domain  of  attraction  of  the  periodic  orbit  of  a  simplified 
three  dimensional  version  of  the  closed  loop  cart  and  pendulum  dynamics  (obtained  by  replacing  the  second 
order  cart  dynamics  with  a  first  order  dynamics).  We  are  also  currently  applying  the  ideas  to  the  nonlinear 
Hoo  control  of  systems  with  periodic  orbits. 
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